An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM (R) and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website. It exhibits good parallel performances (up to similar to 90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. (C) 2014 Elsevier B.V. All rights reserved.

Socio-environmental impact of idol immersion: dimensions and directions

India possesses a diverse and rich cultural heritage and is renowned as a 'land of festivals'. These festivals attract massive community involvement paving way to new materials such as 'Plaster of Paris' being used for 'modernizing' the representation of idols with very little thought given to the issues of toxicity and environmental impacts. Another dimension to the whole issue is the plight of the artisans and the workers involved in the trade. Owing to the unorganized nature of the industry there is minimal or no guidelines pertaining-to the worker safety and health risks of the people involved. This paper attempts to address the complexities of the inherent hazards as a consequence of these socio-environmental issues and trace the scientific rationale in addressing them in a practical and pragmatic way.

Modelling heterogeneity of concrete using 2D lattice network for concrete fracture and comparison with AE study

In this paper, numerical modelling of fracture in concrete using two-dimensional lattice model is presented and also a few issues related to lattice modelling technique applicable to concrete fracture are reviewed. A comparison is made with acoustic emission (AE) events with the number of fractured elements. To implement the heterogeneity of the plain concrete, two methods namely, by generating grain structure of the concrete using Fuller's distribution and the concrete material properties are randomly distributed following Gaussian distribution are used. In the first method, the modelling of the concrete at meso level is carried out following the existing methods available in literature. The shape of the aggregates present in the concrete are assumed as perfect spheres and shape of the same in two-dimensional lattice network is circular. A three-point bend (TPB) specimen is tested in the experiment under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/sec and the fracture process in the same TPB specimen is modelled using regular triangular 2D lattice network. Load versus crack mouth opening isplacement (CMOD) plots thus obtained by using both the methods are compared with experimental results. It was observed that the number of fractured elements increases near the peak load and beyond the peak load. That is once the crack starts to propagate. AE hits also increase rapidly beyond the peak load. It is compulsory here to mention that although the lattice modelling of concrete fracture used in this present study is very similar to those already available in literature, the present work brings out certain finer details which are not available explicitly in the earlier works.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.

Modelling of symmetric laminates under extension

A higher-order theory of laminated composites under in-plane loads is developed. The displacement field is expanded in terms of the thickness co-ordinate, satisfying the zero shear stress condition at the surfaces of the laminate. Using the principle of virtual displacement, the governing equations and boundary conditions are established. Numerical results for interlaminar stresses arising in the case of symmetric laminates under uniform extension have been obtained and are compared with similar results available in the literature.

Modelling the chemical weathering fluxes at the watershed scale in the Tropics (Mule Hole, South India): Relative contribution of the smectite/kaolinite assemblage versus primary minerals

We investigate the chemical weathering processes and fluxes in a small experimental watershed (SEW) through a modelling approach. The study site is the Mule Hole SEW developed on a gneissic basement located in the climatic gradient of the Western Ghats, South India. The model couples a lumped hydrological model simulating the water budget at the watershed scale to the WITCH model estimating the dissolution/precipitation rates of minerals using laboratory kinetic laws. Forcing functions and parameters of the simulation are defined by the field data. The coupled model is calibrated with stream and groundwater compositions through the testing of a large range of smectite solubility and abundance in the soil horizons. We found that, despite the low abundance of smectite in the dominant soil type of the watershed (4 vol.%), their net dissolution provides 75% of the export of dissolved silica, while primary silicate mineral dissolution releases only 15% of this flux. Overall, smectites (modelled as montmorillonites) are not stable under the present day climatic conditions. Furthermore, the dissolution of trace carbonates in the saprolitic horizon provides 50% of the calcium export at the watershed scale. Modelling results show the contrasted behavior of the two main soil types of the watershed: red soils (88% of the surface) are provider of calcium, while black soils (smectite-rich and characterized by a lower drainage) consumes calcium through overall carbonate precipitation. Our model results stress the key role played by minor/accessory minerals and by the thermodynamic properties of smectite minerals, and by the drainage of the weathering profiles on the weathering budget of a tropical watershed. (C) 2010 Elsevier B.V. All rights reserved.

Modelling, Design and Analysis of Low Frequency Platform for Attenuating Micro-Vibration in Spacecraft

One of the most important factors that affect the pointing of precision payloads and devices in space platforms is the vibration generated due to static and dynamic unbalanced forces of rotary equipments placed in the neighborhood of payload. Generally, such disturbances are of low amplitude, less than 1 kHz, and are termed as ‘micro-vibrations’. Due to low damping in the space structure, these vibrations have long decay time and they degrade the performance of payload. This paper addresses the design, modeling and analysis of a low frequency space frame platform for passive and active attenuation of micro-vibrations. This flexible platform has been designed to act as a mount for devices like reaction wheels, and consists of four folded continuous beams arranged in three dimensions. Frequency and response analysis have been carried out by varying the number of folds, and thickness of vertical beam. Results show that lower frequencies can be achieved by increasing the number of folds and by decreasing the thickness of the blade. In addition, active vibration control is studied by incorporating piezoelectric actuators and sensors in the dynamic model. It is shown using simulation that a control strategy using optimal control is effective for vibration suppression under a wide variety of loading conditions.

Modelling of Fe2 + oxidation by Thiobacillus ferrooxidans

The kinetics of oxidation of aqueous acidic ferrous sulphate by Thiobacillus ferrooxidans has been studied in a batch reactor. The contribution of cell wall envelopes to the oxidation rate has been shown to be negligible. A model which accounts for the oxidation of Fe2 +, death of bacteria due to Fe3 + poisoning, existence of an optimal pH and precipitation of Fe3 + has been proposed. The model is able to predict the concentration of Fe2 + and pH quite satisfactorily. The predictions of Fe3 + are not so accurate because of simplifying assumptions made about its precipitation.

Modelling of a batch sonochemical reactor

Ultrasonication of aqueous KI solution is known to yield I2 due to reaction of iodide ions with hydroxyl radicals, which in turn are generated due to cavitation. Based on this conceptual framework, a model has been developed to predict the rate of iodine formation for KI solutions of various concentrations under different gas atmospheres. The model follows the growth and collapse of a gas—vapour cavity using the Rayleigh—Plesset bubble dynamics equation. The bubble is assumed to behave isothermally during its growth phase and a part of the collapse phase. Thereafter it is assumed to collapse adiabatically, yielding high temperatures and pressures. Thermodynamic equilibrium is assumed in the bubble at the end of collapse phase. The contents of the bubble are assumed to mix with the liquid, and the reactor contents are assumed to be well stirred. The model has been verified by conducting experiments with KI solutions of different concentrations and using different gas atmospheres. The model not only explains these results but also the existence of a maximum when Ar---O2 mixtures of different compositions are employed.

Modelling of a batch sonochemical reactor

Ultrasonication of aqueous KI solution is known to yield I2 due to reaction of iodide ions with hydroxyl radicals, which in turn are generated due to cavitation. Based on this conceptual framework, a model has been developed to predict the rate of iodine formation for KI solutions of various concentrations under different gas atmospheres. The model follows the growth and collapse of a gas-vapour cavity using the Rayleigh-Plesset bubble dynamics equation. The bubble is assumed to behave isothermally during its growth phase and a part of the collapse phase. Thereafter it is assumed to collapse adiabatically, yielding high temperatures and pressures. Thermodynamic equilibrium is assumed in the bubble at the end of collapse phase. The contents of the bubble are assumed to mix with the liquid, and the reactor contents are assumed to be well stirred. The model has been verified by conducting experiments with KI solutions of different concentrations and using different gas atmospheres. The model not only explains these results but also the existence of a maximum when Ar-O2 mixtures of different compositions are employed.

Modelling plane mixing layers using vortex points and sheets

We present here a critical assessment of two vortex approaches (both two-dimensional) to the modelling of turbulent mixing layers. In the first approach the flow is represented by point vortices, and in the second it is simulated as the evolution of a continuous vortex sheet composed of short linear elements or ''panels''. The comparison is based on fresh simulations using approximately the same number of elements in either model, paying due attention in both to the boundary conditions far downstream as well as those on the splitter plate from which the mixing layer issues. The comparisons show that, while both models satisfy the well-known invariants of vortex dynamics approximately to the same accuracy, the vortex panel model, although ultimately not convergent, leads to smoother roll-up and values of stresses and moments that are in closer agreement with the experiment, and has a higher computational efficiency for a given degree of convergence on moments. The point vortex model, while faster for a given number of elements, produces an unsatisfactory roll-up which (for the number of elements used) is rendered worse by the incorporation of the Van der Vooren correction for sheet curvature.

Modelling transport processes during microwave heating: A review

The review is concerned with models that analyze transport:processes that occur during microwave heating. Early models on microwave. heating used Lambert's law to describe the microwave power absorption. Over the last decade, models for transport processes have been developed with the microwave power derived from Maxwell's equations. Those models, primarily based on plane waves, have been used for analyzing microwave heating of solids, liquids, emulsions, microwave thawing and drying. The models illustrate phenomena such a resonances, hot spots, edge and runaway heating. The literature on microwave sintering, susceptor heating and microwave assisted synthesis is largely experimental in nature and only key issues are highlighted. To fully appreciate the models for microwave heating, a section on the theory of electromagnetic wave propagation is included, where expressions for the electric field in dielectric slabs and cylinders are presented.

Modelling of a patch repair to a thin cracked sheet

Integrity enhancement of damaged or design deficient structures through repairs is attracting considerable engineering attention. Bonded composite patch repairs to cracked metallic sheets offer various advantages over riveted doubler type, particularly for airframe applications. This paper first reviews the R&D activity in the area of structural repairs. It then approaches the problem of a composite patch repair to a cracked aluminium sheet with different finite element modelling strategies and compares their outcome. The efficient finite element modelling approach thus established is used to study the effect of patch material, patch size, patch symmetry and adhesive thickness on repair performance as the crack grows in the repair configuration. (C) 1999 Elsevier Science Ltd. All rights reserved.